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ABSTRACT 

We measure the z = i?-band optical luminosity function (LF) for galaxies selected 
in a blind Hi survey. The total LF of the Hi selected sample is flat, with Schechter 
parameters M* = — IQ.SS^Qgj + 51og^ioo mag and a — — I.OS^q^^, in good agree- 
ment with LFs of optically selected late-type galaxies. Bivariate distribution functions 
of several galaxy parameters show that the H i density in the local Universe is more 
widely spread over galaxies of different size, central surface brightness, and luminosity 
than is the optical luminosity density. The number density of very low surface bright- 
ness (> 24.0 magarcsec"^) gas-rich galaxies is considerably lower than that found in 
optical surveys designed to detect dim galaxies. This suggests that only a part of the 
population of LSB galaxies is gas rich and that the rest must be gas poor. However, 
we show that this gas-poor population must be cosmologically insignificant in baryon 
content. The contribution of gas-rich LSB galaxies (> 23.0 magarcsec"^) to the local 
cosmological gas and luminosity density is modest (I8I5 ^^"^ ^^2 P^'" '^^^^ respec- 
tively); their contribution to fimatter is not well-determined, but probably < 11 per 
cent. These values are in excellent agreement with the low redshift results from the 
Hubble Deep Field. 
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1 INTRODUCTION 

Understanding galaxy evolution requires well-determined lo- 
cal benchmarks. One of the most fundamental of these is the 
field galaxy luminosity function, the shape of which should 
be predicted by any reliable galaxy formation theory. In 
principle, the shape of the luminosity function is related to 
the power spectrum of primordial density fluctuations and 
complex processes such as gas cooling, star formation and 
feedback to the interstellar medium as well as the behaviour 
of dark matter as it undergoes gravitational collapse and 
merging in galaxy halos (see e.g.. Cole et al. 2000 for a recent 
review). Reference points in the local Universe will help in 
developing a full understanding of these processes. Another 
motive for determination of the local luminosity function, 
is the problem of the faint blue galaxies. The normalisation 
of the z = luminosity function seems to be too low to 
be reconciled with no-evolution predictions based on inter- 
mediate redshift (z ~ 1) surveys (Ellis 1997; Broadhurst, 
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Ellis & Shanks 1988; Koo & Kron 1992), but to quantify 
this problem a reliable measurement of the faint end slope 
is essential. 

The last few years have seen a proliferation of published 
luminosity functions from optical redshift surveys of the lo- 
cal (z < 0.2) Universe (see e.g., Zucca et al. 1997; RatclifFe 
et al. 1998; Folkes et al. 1999; and Blanton et al. (2001) for 
some recent examples). These surveys systematically pro- 
duce samples of lO** galaxies and are able to determine 
the luminosity function down to absolute magnitude limits 
of Mb = —14. However, considerable uncertainty remains 
about the exact shape and normalisation of the luminosity 
function. Especially the faint-end slope for the dwarf galax- 
ies {AIb > —18) is practically unconstrained (see discussion 
in Driver & Phillipps 1996). 

A potential cause of the uncertainty in low z galaxy 
counts is the surface brightness selection effect (Disney 1976; 
Disney & PhiUips 1987). Sprayberry et al. (1997) specifically 
searched for the low surface brightness (LSB) galaxies in 
the APM survey (Impey et al. 1996), and concluded that 
including LSB galaxies in the low z census steepens the field 
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luminosity function, but still does not close the gap between 
number counts at moderate redshift and z = 0. 

The discussion on LSB galaxies ties in directly with an- 
other important benchmark at 2: = 0: the distribution func- 
tion of optical surface brightnesses. Based on photometry 
of 36 nearby spiral and SO galaxies, Freeman (1970) con- 
cluded that ~ 80 per cent have a B band central surface 
brightness /is(0) in the range 21.65 ± 0.30 magarcsec"^. 
The eight deviant galaxies consisted of one dwarf irregular 
LSB galaxy {fisiO) = 23.7 magarcsec"'^) and seven brighter 
galaxies of various morphological type. The majority opin- 
ion at the present moment seems to be that the distribution 
function is flat (McGaugh 1996; de Jong 1996; Dalcanton et 
al. 1997b; O'Neil & Bothun 2000; de Jong & Lacey 2000; 
Blanton et al. 2001; Cross et al. 2001), although Sprayberry 
et al. (1996) found a distribution function that peaks at 
~ 22 magarcsec"^. TuUy & Verheijen (1997) have a dis- 
senting view and present evidence for bimodality in the dis- 
tribution of near-infrared surface brightnesses in the Ursa 
Major Cluster. This view has been contested by Bell & de 
Blok (2000) who claim that the data set is insufficient to 
establish the presence of a bimodal surface brightness dis- 
tribution. 

New insight in both issues can be obtained by selecting 
galaxies via a method that is free from optical selection ef- 
fects. In this paper we measure for the first time the optical 
luminosity function and surface brightness distribution func- 
tion of H I selected galaxies. This sample is the result of the 
Arecibo Hi Strip Survey, a blind strip survey in the 21cm 
line. We stress that this sample is small (60 members) com- 
pared to those produced by modern redshift surveys, and 
large statistical errors are therefore unavoidable. This work 
should be regarded as the first step toward measuring these 
functions for H I selected galaxy samples. Much larger galaxy 
samples will be available in the near future (e.g. HIPASS, 
Staveley-Smith et al. 1996), and the measurements of op- 
tical luminosity functions, surface brightness functions and 
bivariate distributions will greatly improve. 

We organise this paper as follows. First, in section g, we 
briefly describe the sample. In section H we present the op- 
tical luminosity function of this H I selected galaxy sample, 
and discuss the distribution of luminosity density and H I 
gas density among different galaxies. In section H the sur- 
face brightness distribution function and the contribution of 
LSB galaxies to the mass density of the local Universe are 
discussed. Bivariate distribution functions of various galaxy 
parameters are presented in section pi Finally, in section M 
we summarise the conclusions. Throughout this paper we 
use Ho = 100 /iiookms" Mpc~ for calculating distance de- 
pendent quantities. 



2 THE DATA 



column density for the central 3 arcmin wide strip (corre- 
sponding to the main beam of the Arecibo Telescope) was 
~ 10^*cm~^ (5(t) per resolution element of 16 kms~^for gas 
filling the telescope beam. This sensitivity is unmatched by 
any other blind Hi survey to date. Low resolution 21cm 
aperture synthesis observations of the AHiSS sample of 66 
galaxies have been obtained with the NRAO Very Large 
Array (VLA). Details of the Arecibo survey and the VLA 
observations are described by Sorar (1994) and Zwaan et al. 
(1997). 

Optical observations were confined to sources at Galac- 
tic latitudes |fe| > 10° to avoid severe Galactic extinction and 
confusion of foreground stars. This reduces the total number 
of accessible sources to 61. The optical data were obtained at 
the Isaac Newton Telescope of the Observatorio del Roque 
de los Muchachos on the island of La Palma, Spain. The 
data collection was spread over four observing runs during 
the period 1995 October through 1997 March. Images were 
recorded at the Prime Focus camera with a thinned Tek- 
tronix 1024^ pixel CCD. The Tektronix CCD has 24 ^im 
pixels, which give an image scale of 0.59" per pixel at prime 
focus. All images were taken through a standard Harris B 
filter. Flatfields were taken in the twilight, and the residual 
background variations after flatfielding are typically < 1 per 
cent of the sky level. The photometric calibration was done 
by observing standard stars at several airmasses each night, 
and is accurate to 0.13 mag. 

Total galaxy magnitudes were determined using aper- 
ture photometry on the reduced images. Correct aperture 
sizes were found using a curve-of-growth algorithm: aper- 
ture photometry was performed at a series of aperture radii, 
increasing in 1 arcsecond steps, until the integrated mag- 
nitudes levelled out at an asymptotic maximum. The first 
radius at which this maximum (brightest) integrated mag- 
nitude was reached was then chosen as the correct aperture 
size. Central surface brightnesses and disk scale lengths were 
determined by fitting exponential disk models to the az- 
imuthically averaged radial surface brightness profiles. The 
centres of the galaxies were usually taken to be the maxi- 
mum of the light distribution. Strong central concentrations 
were excluded from the exponential fits. The data for galaxy 
A44 turned out to be not usable, which leaves us with a total 
number of 60 galaxies for our analysis. 

We have chosen to apply the internal extinction correc- 
tion proposed by TuUy et al. (1998) which is a function of ab- 
solute magnitude. The extinction correction can be written 
as A^ = 7log(a/fe), where 7 — —0.35(15.1 -I- Mg'^). Since for 
our data set the central surface brightness is well correlated 
with absolute luminosity, the extinction correction implies a 
low correction for LSB galaxies and a higher correction for 
HSB galaxies. Galactic extinction corrections have been ap- 
plied using the reddening maps of Burstein & Heiles (1982) 
and assuming that A{B) =4.1 E{B - V). 



The sample of galaxies used here to measure distribution 
functions of optical luminosity and central surface brightness 
is selected in the 21cm line, and it is therefore free from selec- 
tion effects related to optical surface brightness. The sample 
is the result of the Arecibo H i Strip Survey (AHiSS), a blind 
extragalactic H I survey consisting of two strips of constant 
declination, together covering approximately 65 square de- 
grees of sky over a depth of cz = 7500 kms~^. The limiting 



3 LUMINOSITY FUNCTIONS 

The luminosity function (LF) of galaxies is defined as the 
number of galaxies per cubic Mpc in a luminosity interval 
dM centred at magnitude M. The interval dM is generally 
taken to be 1 mag. The most used parameterisation of the 
luminosity function is the Schechter (1976) function defined 
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by 

<^(M)dM = 0.4 In 10 0* [10 



exp[-lO 



0.4(M*-M)il+c« 



0.4(M*-M)i 



dM, 



(1) 



where a is the faint-end slope, (jj* is the normalisation factor 
and M* is the characteristic absolute magnitude that defines 
the boundary between the exponential and power-law part. 



3.1 Methods 

Many different galaxy luminosity function estimators can 
be found in the literature. In Zwaan et al. (1997) we discuss 
different luminosity function estimators and conclude that 
the E(l/Knax) method is the preferred way to determine 
mass functions and luminosity functions for our sample. For 
this sample, we demonstrated that the determination of the 
Hi mass function with the E(l/Vmax) method is not very 
sensitive to density fluctuations due to large scale structure. 
The E(l/V'max) method consists of summing the reciprocals 
of the volumes corresponding to the maximum distances at 
which galaxies could be seen and still remain within the 
sample. Summing these values per bin in H I mass or abso- 
lute magnitude immediately gives the binned H I mass func- 
tion or optical luminosity function. The advantages of the 
E(l/ymax) method are that it is automatically normalised 
and non-parametric; it recovers the amplitude and the shape 
of the luminosity function simultaneously, without using 
the Schechter function as an assumption about the intrinsic 
shape. An overview of the different galaxy luminosity func- 
tion estimators is given by Willmer (1997), who tests the 
validity of different methods by means of Monte-Carlo sim- 
ulations. Careful examination of his tables shows that the 
E(l/Knax) method (with binning in magnitudes) recovers 
the input luminosity function satisfactorily, and equally well 
as the more conventional parameterised maximum likelihood 
method (Sandage, Tammann & Yahil 1979) or the Stepwise 
Maximum Likelihood Method (SWLM, Efstathiou, Ellis & 
Peterson 1988). Supported by this, we choose to apply the 
E(l/yinax) method to evaluate the optical luminosity func- 
tion. The details of the determination of the values of Vmax 
are described in Zwaan et al. (1997) and will not be repeated 
here. Note that the values of Knax are derived from the origi- 
nal Arecibo H I survey parameters and are independent from 
the optical data. 



3.2 Results 

The resulting luminosity function (^(Ms) is shown in Fig. hi 
as solid dots with la errorbars. The data are binned per 
1.5 mag in order to obtain a reasonable number of galaxies 
per bin, but scaled in such a way that (p represents num- 
ber densities per magnitude bins. Furthermore, the data are 
multiplied by a factor 66/60 to account for the galaxies for 
which no optical information is available. To enable direct 
comparison with published luminosity functions, we choose 
to use the absolute magnitudes uncorrected for opacity ef- 
fects in the galactic disk (see Leroy & Portilla 1998 for a 
discussion on the influence of optical depth effects on the 
shape of the luminosity function). The line indicates the 
best fit Schechter function which is determined by minimis- 
ing x^ for the expected number of galaxies per bin. The 
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Figure 1. Luminosity function for Hi selected galaxies. The 
points were determined using the Vmax method, the errorbars 
are Ic uncertainties from Poisson statistics. The line is the 
best fit Schechter function with parameters: a = — 1.03_Qj^g, 
M* = -19.38lo g2 + 51og/iioo mag and </>* = (1.15 ± 0.40) x 



10" 



noo 



Mpc ^. The inset shows the la and 2a joint two- 



parameter confidence levels for a and M* 



uncertainties in the best fit are indicated in the inset that 
shows the la and 2a error contours of the x^ fit for a 
and M* fitted jointly. As is usually the case in these fits, 
the parameters a and M* are strongly correlated in the 
sense that steeper faint end slopes imply brighter values 
of Af*. The best fit Schechter parameters are found to be 



a = -1.03 



-1-0.25 



M* = -19.38: 



-I- 5 log hioo mag and 



<;/>* = (1.15 ±0.40) X lO"-^ /ij'ogMpc"'', where the quoted er- 
rors are la one-parameter uncertainties. The uncertainties 
given here are solely the result of counting statistics, and 
therefore may understate the true uncertainties. Measure- 
ment errors in the parameters that define Vmax and nrea- 
surement errors in Mb also contribute to the uncertainties, 
but these are relatively small compared to the Poisson errors 
for this small sample. 

The parameterisation in the form of a Schechter func- 
tion is a satisfactory representation of luminosity function 
of the AHiSS galaxies. However, due to the small number 
of galaxies in the low luminosity bins, the value of the faint 
end slope a is poorly constrained. Especially for magnitudes 
fainter than Mb = —14, the slope of the LF is almost un- 
constrained. There is no need for a modification of the LF, 
such as the Schechter function plus a power law, proposed 
by Sprayberry et al. (1997) for his sample of LSB galaxies, 
although our present sample does not rule out this extra 
component. We note that our measured LF parameters are 
in excellent agreement with those from a preliminary anal- 
ysis of the HIPASS survey (Marquarding 2000). 
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3.3 Comparison with optical determinations of 
the LF 

It is interesting to compare tiie luminosity function for H I 
selected galaxies to luminosity functions of optically selected 
galaxies. Recently, there has been much interest in steep 
faint-end slopes of the luminosity function, and the galaxies 
responsible for this steep part 1) are found to be of late 
morphological type (e.g., Marzke et al. 1998), 2) show strong 
emission lines indicative of active star formation (e.g., Zucca 
et al. 1997), and 3) have blue colours (Lin et al. 1999). These 
are the types of galaxies that are expected to contain high 
fractions of H I, and therefore should be represented in the 
AHiSS sample. 

A vast number of luminosity functions based on opti- 
cal redshift redshift surveys, is available in the literature. 
All these surveys contain typically a few thousand galax- 
ies. When making a comparison with our luminosity func- 
tion for Hi selected galaxies, we will concentrate on those 
studies which have made a specific distinction between late 
and early type galaxies, or star forming and quiescent galax- 
ies. We consider: the Stromlo-APM redshift survey (APM, 
Loveday et al. 1992), the Center for Astrophysics redshift 
survey (CfA, Marzke et al. 1994), the ESQ Slice Project 
(ESP, Zucca at al. 1997) , the Las Campanas Redshift Survey 
(LCRS, Lin et al. 1996), the Autofib Redshift Survey (ARS, 
Heyl et al. 1997), the Second Southern Sky Redshift Sur- 
vey (SSRS2, Marzke et al. 1998), the CNOC Field Galaxy 
Redshift Survey (CNOC2, Lin et al. 1999) and preliminary 
results from the 2dF survey (Folkes et al. 1999). 

Table hi summarises the Schechter parameters of 'late 
type' luminosity functions of these surveys. In Fig. H these 
functions are represented, together with our measured points 
and the best fit Schechter function. The functions have been 
transformed to the B filter using the conversions Mb ~ 
Mz - 0.45 for the CfA, 6 - r = 1.1 for the LCRS, and 
Mb — Mbj — 0.24, and all functions are recalculated for 
Ho = 100/iiookms"^Mpc"\ 

The direct comparison of these luminosity functions is 
rather naive for a number of reasons. Firstly, different optical 
wave bands have been used in the selection of these galax- 
ies. This effect can be corrected for by applying a magnitude 
correction, but this is most certainly a oversimplification of 
the problem. The use of different wave bands does not only 
have an influence on the luminosity of the selected galax- 
ies but surely also on the morphological classifications. Sec- 
ondly, the separation between late and early type galaxies 
has been made in different ways for each sample. In the CfA 
sample a detailed separation between morphological types 
has been made on the basis of the galaxies' appearances on 
Palomar Sky Survey. Also the SSRS2 and the APM samples 
have been classified by visual inspection. The CNOC2 data 
is split into different populations using colour information 
of the galaxies. For the 2dF and Autofib surveys spectral 
information has been used to make the classifications. The 
selection criteria for the ESP and the LCRS samples has 
been the occurrence of emission lines in the spectra. In the 
LCRS sample a distinction has been made on basis of the 
criterion of [On] 3727 Wx > 5 A, in the EPS sample the 
selection was simply based on the detection the [On] line. 

With these restrictions in mind, we can compare the 
different luminosity functions for optically selected galax- 




Figure 2. Luminosity functions for late type galaxies. The points 
are the same as in Fig. hi. The lines show the luminosity functions 
from several recent redshift surveys. The details are given in Ta- 
ble hi. Some of these are the summations of several luminosity 
functions for different types. 



ies with Hi selected galaxies. What is particularly striking 
is that the values of the faint-end slope span a wide range 
from -0.80 for the APM survey to ~ -1.50 for the ESP and 
2dF surveys. Even for surveys that use comparable methods 
for classifying their different galaxy population, the differ- 
ences in faint-end slope can be large. Evidently, the shape 
of the luminosity distribution of late type galaxies is still 
ill-constrained. On the other hand, the normalisation and 
the value that defines the knee are quite similar for all sur- 
veys; all functions cross approximately the same point at 
Mb ~ —19-1- 51ogftioo. The luminosity function for the 
AHiSS falls in between those of the optical samples. We 
therefore conclude that our estimate of the luminosity func- 
tion is in good agreement with that of optically selected 
samples. Furthermore, there is no new population selected 
by H I surveys that adds significantly to the galaxy popula- 
tions identified through optical surveys. 



3.4 Luminosity density of gas-rich galaxies 

A more fundamental parameter is the luminosity density, 
the integrated light from the whole population of galaxies. 
As is discussed by Lilly et al. (1996), this parameter is in 
principle less dependent on the details of galaxy evolution 
than the luminosity function. The integral luminosity den- 
sity of late type galaxies can be determined by integrating 
the Schechter luminosity function weighted by luminosity, 
which gives jb ~ (j>* L*bT{2 + a), where F is the Euler 
gamma function. The values of jb for late-type galaxies as 
determined by the different optical surveys is given in the 
last column of Table |l| It is remarkable that all values of 
JB are within ~ 1.5(t from the value determined from the 
AHlSS. A notable exception is the 2dF survey that finds a 
value 60 per cent higher than the mean of the other surveys. 
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Table 1. Comparison of Luminosity Functions for Late Type Galaxies 



Sample 


Selection 


a 


M* -Sloghioo 


0*- 


Pl"^ 


AHiSS 


H I selected 


-1.03 


-19.38 


11.5 


10.3 ±2.0 


APM (Loveday et al. 1992) 


Sp/Irr 


-0.80 


-19.16 


10 


6.7 


CfA (Marzke ct al. 1994) 


Sa-Sb 


-0.58 


-18.93<= 


8.7 


8.9 




Sc-Sd 


-0.96 


-19.02<= 


4.4 






Sm-Im 


-1.87 


-19.00<= 


0.6 




ESP (Zucca et al. 1997) 


Emission lines 


-1.40 


-19.23 


10 


11.4 


LCRS (Lin et al. 1997) 


3727 W^A > 5 A 


-0.90 


-18.93'^ 


13 


7.7 


Autofib (Heyl et al. 1997) 


Sab 


-0.99 


-19.76 


2.19 


8.4 




Sbc 


-1.25 


-19.16 


2.80 






Scd 


-1.37 


-18.96 


3.01 






Sdm 


-1.36 


-18.76 


0.50 




SSRS2 (Marzke et al. 1998) 


Spirals 


-1.11 


-19.43 


8.0 


9.2 




Irr/Pec 


-1.81 


-19.78 


0.2 




CNOC2 (Lin et al. 1999)<= 


Intermediate type 


-0.53 


-18.97 


9.0 


11.1 




Late type 


-1.23 


-19.07 


7.2 




2dF (Folkes et al. 1999) 


Sab 


-0.86 


-19.44 


3.9 


14.0 




Sbc 


-0.99 


-19.14 


5.3 






Scd 


-1.21 


-18.76 


6.5 






Sdm-Im 


-1.73 


-18.78 


2.1 




'^ Units are 10"^ /ifpg Mpc'^ 












'^ Units are lO'^ hioo Lg Mpc" 


-3 










" B - Mz = -0.21 












d B-_R= 1.1 












° Values extrapolated to 2 = 














Folkes et al. (1999) note that the 2dF results are prelimi- 
nary, and that corrections for completeness, clustering, and 
Malmquist-bias have not been applied yet. It remains to be 
seen whether the final 2dF results will remain in excess of 
the AHiSS estimate of jb- A preliminary result from the 
SDSS (Blanton et al. 2001) has a substantially higher op- 
tical luminosity density than other recent surveys, but this 
increase appears to arise in a their photometric evaluation 
of each galaxy's luminosity rather than an increase in the 
number density of objects. 

The mean value of jb of all optically selected late- 
type galaxy samples is 9.7 x 10^ /iioo L© Mpc~^, while that 
for the AHiSS sample is (10.3 ± 2.0) x lO'^ /iioo L| Mpc'^ 



This latter value translates to js 



W 



'/irj'oWHz-'Mpc- 



(3.4 ± 0.7) 



, using the conversion of Lilly et al. 



(1996). This is approximately 50 per cent of the integral lu- 
minosity density of the local Universe as measured by the 
most recent optical redshift surveys using isophotal magni- 
tudes (Folkes et al. 1999, Blanton et al. 2001). Blanton et 
al. (2001) find that jb increases significantly if extrapolated 
magnitudes are used. 



3.5 Luminosity and Hi mass distributions for 
different morphological types 

A more detailed view of the relative importance of different 
morphological types to the Hi and luminosity density can 
be made be transforming luminosity functions into H I mass 
functions, assuming correlations between Hi mass and op- 
tical luminosity. Rao & Briggs (1993) used this method to 
determine the Hi mass function and fiHi based on at that 
time available luminosity functions. They showed that by 



adopting the relation log Mhi = a — bMs between H I mass 
and optical luminosity, the H I mass function can be written 



e(MHi)d(MHi) = ^^*(|||i)'"+^'^-x 



(2) 



where log Mhi = a — &M^ , and a and 0* are the Schechter 
parameters of the luminosity functions. 

Here we update the calculations by Rao & Briggs (1993) 
with more recent luminosity function parameters, and test 
if the results are in agreement with our measurements. For 
completeness, we present all possible ways of plotting the 
number density, the H I density, and the luminosity density 
as a function of absolute magnitude and H I mass. We adopt 
the Marzke et al. (1998) luminosity functions for different 
morphological types, and we fit linear regression lines to 
Mhi vs. Mb taken from the Nearby Galaxy Catalog (Tully 
1988) to find the values of a and b. 

Fig. H shows the results. The solid grey lines are for all 
galaxy types, the black dashed lines are for spirals, the grey 
dashed lines for E and SO types, and the grey dotted line for 
Irr/Pec types. The thin parts of each line are extrapolations 
beyond the confidence levels set by Marzke et al. (1998). The 
solid points plus errorbars are different representations of 
the AHiSS data. The results are basically the same as what 
Rao & Briggs (1993) found. We show here that the optical 
luminosity functions, combined with conversion factors from 
Mb to Mhi , give excellent fits to our data. The H i density 
distribution matches the converted luminosity functions for 
all galaxy types summed, and the luminosity distribution 
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Figure 3. Luminosity and H I density functions for different mor- 
phological types. The points are from the AHiSS and are the re- 
sult of the l/Vmax method. The lines are converted functions for 
different morphological types, using Marzke et al. (1998) luminos- 
ity functions, and Hi to Mb relations from the Nearby Galaxy 
Catalog (TuUy 1988). The solid grey lines are for all galaxy types, 
the black dashed lines for spirals, the grey dashed lines for E and 
SO, and the grey dotted lines for Irr/Pec. The thin parts of each 
line are extrapolations beyond the confidence levels set by Marzke 
et al. (1998). Upper left panel: Luminosity function. Upper right 
panel: Hi mass function. Middle left panel: Luminosity density 
as a function of Mg. Middle right panel: Luminosity density as 
a function of Mjji. Lower left panel: Hi density as a function of 
Alg. Lower right panel: Hi density as a function of Mjji. 



is fitted satisfactorily witli spiral and irregular population. 
It is no surprise that the luminosity density from ellipticals 
exceeds the measured values from the AHiSS, since these 
objects are not selected by Hi surveys. This is especially 
clear in the middle right panel. 

The integral Hi density can be determined from a op- 
tical luminosity functions via 



PHI 



$(L) Mm dL = (j)* 10"-''" r{l + a + 2.5 &).(3) 



If we apply this to our adopted luminosity functions we 
find that spirals make up 62 per cent of the Hi gas den- 
sity, Irr and Pec types contribute 35 per cent, and E and 
SO types only 3 per cent. Natarajan & Pettini (1997) apply 
this same method to measurements of the luminosity func- 
tion at higher redshift in order to chart the evolution of the 



cosmic gas content between z = 1 and z — 0. The viability 
of this result is unclear since the amount of evolution of the 
Mhi/L ratio of galaxies is presently unknown. Future deep 
Hi surveys at redshifts z > are required to constrain the 
Mhi/L evolution. 



4 CONTRIBUTION OF LSB GALAXIES TO 
THE COSMIC MASS BUDGET 

More than two decades after the seminal paper by Disney 
(1976) who defined the potential selection effects against 
LSB galaxies, the debate on the cosmological significance of 
LSB galaxies is still open. The AHiSS sample, which is not 
biased by the sky background, makes a valuable contribu- 
tion to the discussion of the cosmological significance of LSB 
galaxies. 

4.1 The surface brightness distribution function 

The bottom panel of Fig. H shows a histogram of the distri- 
bution of _B-band surface brightnesses for the AHiSS galax- 
ies. These surface brightnesses are corrected for dust extinc- 
tion following the formalism described by TuUy et al. (1998). 
The unshaded histogram shows the distribution for the full 
set of AHiSS galaxies, and the grey histograms show the 
distribution for the subset of AHiSS galaxies with inclina- 
tions i < 75° for which the corrections to face-on values are 
modest. The galaxies with high inclinations {i > 75°) do 
not appear from their optical images to be very low surface 
brightness: they often exhibit bright central condensations 
and strong dust lanes, both features not normally found in 
extreme LSB galaxies (see McGaugh, Schombert, & Bothun 
1995). It therefore seems possible that the true face-on sur- 
face brightness of these disks is brighter than those given by 
either the TuUy et al. (1998) prescription or the assumption 
that the disks are fully transparent. 

The top panel of Fig. Ml shows the volume-corrected sur- 
face brightness distribution function of AHiSS galaxies. This 
function is determined by summing values of 1/Vmax per 1 
mag bins of surface brightness. The errorbars indicate 68 
per cent confidence levels and are determined from 100 boot- 
strap re-sample realizations of the data. The hollow symbols 
show the distribution for the complete set of AHiSS galaxies, 
and the solid symbols are limited to those AHiSS galaxies 
with i < 75°. Note that the distribution function resembles 
the one found by Sprayberry et al. (1996) based on the APM 
survey. 

4.2 A cutoff in surface brightness? 

As Fig. W shows, the AHiSS detected no galaxies with reli- 
ably determined face-on central surface brightnesses fainter 
than fj,g^ = 24 magarcsec"^. Even among the highly in- 
clined galaxies with large (and possibly unreliable) correc- 
tions to face-on values, there are no galaxies with /is ^ 
25 magarcsec"^. The statistical significance of this results 
depends on the assumptions we make about the detectability 
of very LSB systems. We can make an estimate by calculat- 
ing the average value of Vmax for the different surface bright- 
ness bins. We find that Knax is mildly correlated with p^'': 
dimmer galaxies can on average be detected over smaller 



© 0000 RAS, MNRAS 000, §-^ 



The Luminosity Function and Surface Brightness Distribution of Hi Selected Galaxies 7 



OB 

CO 



Si 




3- 




on 




() 






-b 




?.0 






ijj 


15 


6 

3 


10 


C 


5 








Dalcanton et ol. (1997) 



' \ ' r 

o All 

• inci < 75° 



f y 




Figure 4. Bottom panel: Distribution of surface brightnesses in 
the AHiSS sample. The grey histogram is for those galaxies with 
i < 75° and is embedded in the histogram for all galaxies. Top 
panel: Volume corrected distribution of surface brightnesses. The 
open and solid symbols have the same meaning as in the bottom 
panel. For clarity, the points are slightly offset horizontally. Er- 
rorbars indicate 68 per cent confidence levels. Arrows denote 95 
per cent confidence upper limits. The space density of optically 
selected LSB galaxies determined by Dalcanton et al. (1997b) is 
indicated by a light grey box, and corresponds to 90 per cent 
confidence levels. 



volumes. This correlation arises because decreasing surface 
brightness correlates with decreasing total Hi mass, and 
the sample selection is based on H I flux. If we extrapolate 
the iig — Vmax correlation to the surface brightness bins 
in which we have no detections, we find that (Knax) would 



max/ 
2 



be 1150/ij~QgMpc for the 24 — 25 magarcsec" bin, and 
950 h^QQ Mpc^ for the 25 — 26 mag arcsec"^ bin. The proba- 
bility pk of finding k objects when the mean is n, is given by 
the Poisson distribution: pk = e~"n''/k\. The mean num- 
ber of detected objects per mag is given by (p(/i)Vmax(M)j 
where ip{fJ,) is the space density of objects as a function 
of surface brightness. Hence, the probability of finding zero 
sources in one bin is po = e~'^'^^^"'"=''''. A 95 per cent 
confidence upper limit to ip{j-L) can now be expressed as 
ifii^j.) — — ln(0.05)/Vmax(A*)- This equation is used for the 
upper limits that are indicated by arrows in Fig. 0. 

This absence of extreme LSB galaxies suggests two 
things. First, the space density of massive, gas-rich, ex- 
tremely LSB disks such as Malin 1 (Bothun et al. 1987) must 
be low, as previously shown by e.g., Briggs (1990, 1997), 
Driver & Cross (2000) and Blanton et al. (2001). Such a 
disk would have been easily detectable to the limit of 7500 
kms~^ of the AHiSS, so they must be intrinsically less com- 
mon on average than 1 per 1000 /i^qq Mpc^ (95 per cent con- 
fidence level, using the Poisson statistics). Second, optical 
surveys for LSB galaxies (Sprayberry et al. 1996; Dalcanton 
et al. 1997b; O'Neil et al. 1997) systematically find galaxies 
at lower surface brightnesses than ^(0) = 24 magarcsec"^. 



so there must be some reason why the present Hi survey 
fails to detect any. 

One possibility is that such galaxies have detectable 
amounts of neutral hydrogen but are extremely rare, so that 
it would not be expected to find one in the AHiSS search 
volume. Apart from a few special objects like Malin 1, this 
seems unlikely because optical surveys find these objects in 
significant numbers despite the relatively small volume lim- 
its imposed by optical surface brightness selection effects 
(McGaugh 1996). Specifically, Dalcanton et al. (1997b) find 
that the number density of galaxies with V-band central 
surface brightnesses in the range 23 < /x < 25 magarcsec"^ 
is 0.081004 ''■lOO Mpc~'^ comparable to the number density 
of normal galaxies. For reference, we have indicated this es- 
timated with a shaded box in Fig. W, where we have adopted 
B ~V = 0.5, a typical value for LSB galaxies (de Blok, van 
der Hulst & Bothun 1995). 

The other possibility is that a significant number of 
galaxies exist in the AHiSS search volume with optical sur- 
face brightnesses /i(0) > 24 magarcsec"'^, but that they do 
not contain enough Hi to be detected by the AHlSS. This 
seems considerably more likely, as there are two ways such 
a population could come to exist: First, these very low den- 
sity systems could have formed a first generation of stars 
and then either lost most of their remaining gas through 
supernova-driven winds (Babul & Rees 1992; Babul & Fer- 
guson 1996) or consumed all their gas in vigorous star for- 
mation and since then faded to become LSB disks (Bell et 
al. 1999). Second, like the outskirts of normal spiral galax- 
ies, LSB disks have low Hi surface densities (de Blok, Mc- 
Gaugh & van der Hulst 1996) , and as such, they are subject 
to ionization by the extragalactic UV background that pro- 
duces the sharp cutoffs seen at the edges of normal spirals, 
for column densities below 10^®'® cm~^ (e.g., Maloney 1993; 
CorbeUi & Salpeter 1994; Dove & ShuU 1994). Thus, much 
of the gas in LSB disks should become ionised, and thus be 
undetectable in 21cm surveys. This is consistent with the 
finding that no AHiSS galaxies were found with average H I 
column densities lower than {Nui) > 10^®'^ cm~'^ (Zwaan et 
al. 1997). The average limiting column density of the AHiSS 
was ~ 10^*cm~^ (5cr) per 16 kms~^. Of course, for many 
galaxies in the AHiSS sample this number is not the min- 
imal detectable Hi column density averaged over the gas 
disk. Not all detected galaxies fill the beam of the Arecibo 
Telescope, and the velocity width of all detections is larger 
than 16 kms"'^. A typical AHiSS galaxy fills 50 per cent 
of the beam and has a velocity width of 160 kms~^. With 
optimal smoothing applied, the minimal detectable column 
density of such a galaxy would be 10^*'* cm"'^, still approxi- 
mately an order of magnitude lower than the reported cut-off 
at lO^^^'cm"^ in Zwaan et al. (1997). Moreover, the detec- 
tion limit of very large galaxies or gas clouds that do fill 
the beam and have a similar velocity width of 160 kms~^ 
would be even lower, approximately lO^^'^cm"^. None of 
these extended, low column density systems were detected. 

At present there is insufficient data to distinguish be- 
tween the two proposed hypotheses. Currently available 
studies of the stellar compositions of LSB galaxies (Mc- 
Gaugh & Bothun 1994; de Blok et al. 1995; BeU et al. 1999 
and 2000) have concluded that gas-rich LSBs form stars 
slowly and continuously and therefore have fairly young stel- 
lar populations. Judging from their colours, the newly iden- 
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tified class of red LSB galaxies (O'Neil, Bothun & Schombert 
2000) are consistent with a scenario in which they are sim- 
ply a fading, passively evolving population (Bell et al. 1999). 
However, O'Neil et al. (2000) report that some of these 
galaxies have high values of Mni/L, but they also note that 
60 per cent of their LSB galaxy sample that was foUowed- 
up with Arecibo is undetected in H I, in sharp contrast to 
the high success rate in earlier LSB samples (e.g., Schnei- 
der et al. 1990) . Moreover, a cross-correlation of the tables 
in O'Neil et al. (1997) and O'Neil et al. (2000) shows that 
the global V — I colours of the undetected galaxies are on 
average 0.3 mag redder than the galaxies in which Hi was 
found. 

Multicolour photometry of galaxies with ^b > 
24.0 magarcsec"^, in combination with deep H-a imaging 
and deep 21cm observations should show whether the low- 
est surface brightness galaxies are consistent with a fading, 
passively evolving population. 



4.3 The LSB contribution to the neutral gas 
density 

To address the problem of the cosmological significance of 
gas-rich LSB galaxies in a meaningful way, an LSB galaxy 
should be well-defined. In the literature, different authors 
adopt different definitions for the critical surface brightness 
that separates galaxies into the 'normal' and LSB classes. 
The critical value ranges from 21.65 mag arcsec"'^ (the 'Free- 
man value') to 23.5 magarcsec^^. In the remainder of this 
paper we define an LSB galaxy as a galaxy with de-projected 
B-band central surface brightness > 23.0 mag arcsec"'^ . This 
limits the LSB galaxies to those that are ~ 4a dimmer than 
the Freeman value. 

The cumulative distribution of Hi density among 
AHiSS galaxies of different surface brightness is presented 
in the top panel of Fig. g. The H I density distribution can 
be fit satisfactorily with a Gaussian distribution. There is no 
fundamental physical motivation for using a Gaussian to pa- 
rameterise the distribution function, but Dalcanton, Spergel 
& Summers (1997a) note that a galaxy formation scenario 
based on a log-normal distribution of the spin parameter 
A, produces a (nearly) Gaussian function of luminosity den- 
sity vs. surface brightness (see also de Jong & Lacey 2000). 
The inset in Fig. h gives la and 2a confidence ellipses for 
the Gaussian fits, the horizontal axis shows the centre of the 
distribution, the vertical axis the dispersion (la). The actual 
fitting was done on the binned data, not on the cumulative 
distribution. 

The Hi mass density of the local Universe is domi- 
nated by galaxies dimmer than the Freeman (1970) value 
of 21.7 magarcsec"'^. The peak of the differential distribu- 
tion is at 22.0 magarcsec"^, the width is 1.0 magarcsec"^. 
Low surface brightness galaxies contribute a minor fraction 
to the Hi density, galaxies fainter than 23.0 magarcsec"^ 
make up ISlg per cent of the Hi mass density in the lo- 
cal Universe (the quoted errors have been determined using 
bootstrap re-sampling and mark the 68 per cent confidence 
levels) . 
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Figure 5. Cumulative distributions of Hi mass density (top 
panel) and luminosity density (bottom panel) as a function of 
central surface brightness for the AHiSS sample. The lines show 
Gaussian fits. The la confidence levels on the two jointly fitted 
parameters fi and c, the mean and the width of the Gaussian, 
are shown in the inset. The Hi density of the local Universe is 
dominated by galaxies with S-band central surface brightness 
of 22.0 magarcsec"'^, the luminosity density is dominated by 
21.2 magarcsec" galaxies for this Hi selected sample. 



4.4 The LSB contribution to the luminosity 
density 

The lower panel of Fig. ra shows the cumulative distribution 
of luminosity density against surface brightnesses. The peak 
of the differential distribution is at 21.2 magarcsec"^, the la 
dispersion is 1.0 magarcsec"^. This implies that while most 
of the H I resides in galaxies dimmer than the Freeman value, 
most of the light in gas-holding galaxies in the local Universe 
is in galaxies 0.5 mag brighter than the Freeman value. The 
contribution of LSB galaxies is insignificant; galaxies with 
fi^' > 23 magarcsec"^ constitute no more than 5 ± 2 per 
cent to the luminosity density. 

We stress that this result only holds for gas-rich LSB 
galaxies. The contribution of gas-free LSB galaxies is un- 
constrained by our survey. Sprayberry et al. (1997) con- 
clude that optically selected LSB galaxies contribute about 
30 per cent to the field galaxy luminosity density, a result 
very consistent with ours, since their definition of an LSB 
galaxy is fis > 22.0 magarcsec"^. De Jong & Lacey (2000) 
find that the luminosity density of optically selected galax- 
ies is dominated by ^i ~ 19.3 magarcsec"^, which compares 
to 21.0 magarcsec"^ in the B-band (using their value of 
B — I — 1.7). They estimate that approximately 4 per cent 
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of the luminosity density is provided by galaxies with ^,g > 
22.75 magarcsec"^. Driver (1999) defines a volume limited 
subsample of 47 galaxies at 0.3 < z < 0.5 from the Hub- 
ble Deep Field and derives that LSB galaxies (mean surface 
brightness within the effective radius > 21.7 magarcsec"^) 
contribute 7 ± 4 per cent to the luminosity density. All these 
results are in good agreement with our estimates. This fact 
implies that if a population of very LSB, gas-free LSB galax- 
ies exists (as was discussed in Section 1.2), their contribution 



to the luminosity density must be negligible. 



4.5 The LSB contribution to fimattcr 

The contribution of LSB galaxies to the total mass budget 
of the local Universe is a longstanding question. The cal- 
culation critically depends on the assumptions one makes 
on the dependence of the dynamical M/L on central sur- 
face brightness. A zeroth order approximation is to assume 
that M/L is equal for all galaxies, independent of central 
surface brightness. This assumption follows naturally from 
the observation that surface brightness is not a parameter in 
the TuUy-Fisher relation (Sprayberry et al. 1995; Zwaan et 
al. 1995; Verheijen 1997). Moreover, Verheijen (1997) shows 
that it is possible to use one model for the dark matter halo 
to fit the rotation curves of three galaxies, all at equal posi- 
tion in the Tully-Fisher relation, but with different surface 
brightness. In the terminology of McGaugh & de Blok (1998) 
this invariant M/L would be the 'same halo hypothesis.' It 
is consistent with the idea that all galaxies of equal luminos- 
ity form in the same mass halo, but the angular momentum 
of an LSB disk is higher, which causes the disk to be less 
centrally concentrated (Dalcanton et al. 1997a). Comb ined 



with the result on the luminosity density from section 4.4 



this assumption leads to the conclusion that pm(LSB) is 5 
per cent of the total pM (i.e., equal to the contribution to 
the luminosity density). 

Van den Bosch & Dalcanton (2000) show that their 
semi-analytical galaxy models are consistent with M/L ex 
J1-1/2 ^2waan et al. 1995), where E is the central surface 
brightness in linear units. M/L ratios are calculated via 
M oc DV^ , where V is the maximum rotational velocity 
and D is a characteristic size of the dark halo, which is as- 
sumed to be directly proportional to the scale length of the 
optical disk. If we, like Driver (1999), adopt this relation for 
M/L we find that the LSB contribution to pM rises to III3 
per cent. If we apply the calculations of M oc DV^ directly 
to our AHlSS data set, we find pm(LSB) — IOI3 P®'' cent. 
Both values are in excellent agreement with the 12 ± 6 per 
cent that Driver finds. 

At present it is unclear what the true dependence of 
the dynamical M/L on optical surface brightness is. Clearly, 
high precision measurements of rotation curves of LSB sys- 
tems are needed (see Swaters, Madore & Trewhella 2000; 
van den Bosch & Swaters 2000; van den Bosch et al. 2000). 
At the moment we adopt as a conservative estimate that 
gas-rich LSB galaxies contribute no more than 11 per cent 
to pM, the dynamical mass contained in galaxies. 



5 BIVARIATE DISTRIBUTIONS 

A more detailed view of the distribution of baryons among 
galaxies of different size and brightness can be obtained by 
calculating bivariate distribution functions. The importance 
of this way of looking at galaxy parameters is stressed by 
van der Kruit (1987, 1989), de Jong (1996), and most re- 
cently by de Jong & Lacey (2000) who study 10^ galaxies 
with types Sb to Sdm. First, the bivariate distribution func- 
tion is a important constraint for galaxy formation theories, 
as any theory should not only produce the integrated lu- 
minosity function (and integrated distribution functions of 
other structural parameters), but also higher dimensional 
distribution functions. Second, bivariate distribution func- 
tions help to clarify the selection effects that influence the 
determination of (e.g.) the luminosity function. 

The aim of the present work is 1) to test whether an 
Hi selected galaxy sample yields qualitatively the same bi- 
variate distribution function, and, 2) extend the bivariate 
distribution functions to the distribution of Hi properties. 
The sample we study here is small, and, as is discussed by de 
Jong & Lacey (1999) and Minchin (1999), at least 500-1000 
galaxies are required to avoid problems with small number 
statistics. We only intend to make qualitative comparisons 
and care should be taken with the interpretation of the re- 
sults. We also stress that the sample that we use here might 
be biased against galaxies with very low values of Mui/L, 
just like optical samples are biased to those with high values 
of Mm/L. 



5.1 Results 

In Fig. |6| we present bivariate distributions of several fun- 
damental parameters. Our aim is to show the distribution 
of H I mass density and luminosity density as a function of 
galaxy luminosity, gas mass, size, and surface brightness. 

The figures are calculated by distributing values of 
Afni/Vmax and LgYKnax over a fine grid with 0.1 dex resolu- 
tion for A/hi and Lg\ 0.1 mag resolution for central surface 
brightness fig\ and 0.05 dex resolution for disk scale length 
h. Next, the images were smoothed with a Gaussian filter, 
of which the FWHM values are indicated by the dashed el- 
lipses in the upper left corners of each panel. The H I density 
distributions are shown as contours, the luminosity density 
distribution as grey-scales. Steps in intensity are in logarith- 
mic intervals of 0.25 dex. 

The first thing to notice is that the general trends are 
the same for the gas density and the luminosity density, but 
the maximum of the gas density is shifted towards less lu- 
minous, lower surface brightness galaxies. The second point 
is that the luminosity density is more strongly concentrated 
towards large, luminous HSB galaxies, whereas the Hi den- 
sity is more widely distributed. 

What is obvious from the top right panel is that both 
the luminosity and the H I mass distribution are strongly de- 
pendent on optical surface brightness in the sense that both 
functions are shifted towards fainter absolute magnitudes 
for lower surface brightness galaxies. This fact was also ob- 
served by de Jong (1996), and the peak of the distribution, 
at Mg'^ = —20-1-5 log /iioo and pg' = 21 magarcsec"^ agrees 
well with his determination. The same trend is observed by 
Cross et al. (2001) who present a bivariate distribution func- 
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Figure 6. Bivariate distribution of luminosity density (grey-scales) and Hi density (contours) for Hi selected galaxies in the (Lg^, 
MHi)-plane {top left), the {Lg , /jg'')-plane {top right), the {h, /ig'')-plane {lower left), and the {Lg , /i)-plane {lower right). Grey-scales 



correspond to {10^-^ ,10^-''^ ,10^'^ , ...) X hioo i© Mpc-^, contours to (lO'^-^S, lO^-S, 10^ 



.) X /iioo Af0 Mpc ^. The densities are per 



decade for L^'' and Mui, per 0.3 dex for h and per mag for lJ.g ■ The dashed ellipses in the upper right corners of each panel indicate 
the FWHM of the Gaussian smoothing filter that has been applied to the data. 



tion based on a a preliminary subsaniple of 5 x IC* galaxies 
from the 2dF survey and by Blanton et al. (2001) who use a 
sample of 10'* galaxies from the SDSS commissioning data. 
A similar effect can be seen in the lower right panel: the 
luminosity and the H I mass distribution are shifted towards 
fainter absolute magnitudes for smaller galaxies. This corre- 
lation has been studied in detail by de Jong & Lacey (2000) , 
who discuss the predictions of hierarchical galaxy formation 
theories and conclude that the observed distribution is in 
qualitative agreement with theory, but the distribution in 
disk size is narrower than predicted. 



6 CONCLUSIONS 

We have presented several volume corrected distribution 
functions for galaxies that have been selected in the H l 21cni 
line. The conclusions are the following: 



1. The luminosity function of the Hi selected galaxies is 
in agreement with other determinations based on late-type, 
or star- forming galaxies. The integral luminosity density of 
gas-rich, late-type, or star-forming galaxies is well deter- 
mined and equals js = (3.4±0.7) x 10*^ h'^^^ W Hz"* Mpc"^. 
This is approximately 50 per cent of the integral luminosity 
density of the local Universe. 

2. The contribution of low surface brightness (LSB) 
galaxies to the integral luminosity density and Hi den- 
sity is modest, 5 and 18 per cent, respectively. This in in 
good agreement with calculation based on optically selected 
galaxies. The fraction of ^matter that resides in LSB galaxies 
is at present not well determined, but probably less than 11 
per cent. 

3. We observe a lower limit to the surface brightness of 
gas-rich galaxies: no galaxies were found with de-projected 
central surface brightness > 24.0 magarcsec"^ in the B- 
band. It will be interesting to test whether this result stands 
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up with future large H I surveys, such as the HIPASS survey 
(Staveley-Smith et aL 1996). 

4. Bivariate distributions of various fundamental galaxy 
parameters show that the H I density in the local Universe is 
more diffusely spread over galaxies with different size, sur- 
face brightness, and luminosity than the luminosity density. 
The luminosity density is concentrated towards bright, large, 
high surface brightness disks. 
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